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ABSTRACT 

Energy spectra of nearly stripped carbon and oxygen nuclei observed 
during several solar particle events indicate a systematic deviation of 
these spectra from a simple power law: the spectra bend down below about 

100 keV per nucleon and the degree of turn-over is highly correlated with 
the size of the flare, as measured by the "event averaged" flux of 
130-220 keV protons. The energy spectra of helium computed for the same 
time periods do not show a similar featiire. We find a large variability 
of the a/CNO ratio from event to event (from 2 to about 20 at Uo keV per 
nucleon), and in all cases examined the carbon and oxygen nuclei are nearly 
fully stripped. These results are interpreted as evidence for storage of 
energetic ions in hot (T^ 'v 1,5 x 10 °K) coronal regions, followed by 
strong adiabatic deceleration. 


Sub.lect headings: Abundances, solar cosmic rays flares, solar 
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I. INTRODUCTION 

I 

The first n^asurements of the ionization states and energy spectra of 
solar particles extending to energies of ''*20 keV per nucleon have been 
reported recently by Scieuabi et al. (197^) • These results indicated that 
in one moderate-sized solar particle event examined (U-8 July 197^) low 
energy csurbon and oxygen ions were nearly fully stripped and that their spectra 
deviated significantly from a simple power law in kinetic energy. In this 
paper, we describe the results of analysis of eight additional events, and 
conclude that the large degree of ionization and turn-over in the spectra are 
consistent features of low energy solar CNO ions. We examine the roles that 
coronal storage and interplanetary propagation might have in producing the 
observed spectra. 

II. INSTRUMENTATION AND DATA ANALYSIS 

The University of ^feryland experiment aboard the earth-orbiting 
spacecraft IMP-8 (Explorer 50) is an electrostatic deflection sensor 
(Turns ^ , 197^). A simplified cross-sectional view of the instrument 

is shown in Figure 1. An ion entering the instrument has its entrance 
trajectory defined by a collimator. The particle then enters one of three 
electrostatic field regions where its path is deflected by an amount 
inversely proportional to its kinetic energy per charge. On exit from the 
deflection cavity, the ion strikes one of an array of eight rectangular 
solid state detectors, each of which covers a well-defined deflection, and 
hence energy per charge, range. Pulse-height analysis of the detector 
signal allows a direct measurement of the particle's total energy, from 
which its charge may be determined. In addition, by using multiple energy 
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thresholds on signals of several detectors, twelve counting rates having 
well-defined energy limits and hence charge responses are generated with 
time resolutions as short as five seconds. Nine of these rates are 
sectored into four directional quadrants in the ecliptic plane* providing 
data on particle anisotropies. Solid-state anticoincidence detectors, 
placed "behind the detector array, and a plastic scintillator anticoincidence 
cup prevent analysis of penetrating particles. The energy ranges and charge 
responses of the various detectors are shown in Table I. 

Background correction of a detector's pulse-height histogram is 
performed by subtracting a normalized histogram obtained from the same 
detector during quiet times, when the counting rates are at their lowest 
values. One of the solid-state detectors (P5) in the detector array is 
covered by a thin metal foil which stops deflected particles. It thus 
serves as a background monitor which, when pulse-height analyzed, 
identifies any changes in the shape of the background with time. Analysis 
of pulse-height histograms of this background detector shows that during 
quiet times all of the other detectors respond primarily to the secondary 
electron background produced in the spacecraft. 

The relative abundance of individual charge states may be determined 
by fitting computed instrument response curves to the observed pulse-height 
distributions. The sensor is capable of fully resolving ionic charge 
states 1 to 4 for elements up to oxygen. In addition, charge states 5 to 
8 are separated sufficiently so that, by means of a sinple iteration 
procedure, abundances of these individual states may also be determined. 
Complete charge state identification for carbon and oxygen down to about 
liO keV per charge is thus obtained. Of course, the individual elements 
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may not be resolved in this majiner since only the ionic charge (as opposed 
to nuclear charge) is determined. In many cases, however, the ions may be 
completely identified (i.e., ionic charge and nuclear charge) by using the 
simultaneous nuclear abundance information provided by the ULET, a thin- 
window proportional coiinter dE/dX ^ E telescope (Hovestadt and Vollmer, 
1971 ) of the Maryland/Max-Planck-Institut experiment. 

Figure 2a shows a typical uncorrected pulse-height spectrum (histogram) 
obtained during the solar flare particle event of 197^ June 22. The 
measured smoothed background distribution is shown as a dashed line, and 
the histogram after background correction is shown in Figure 2b. This 
particular histogram was obtained from detector P3, which responds to 
ions in the energy range 130 to 220 keV per charge. The triple-peaked 
shape of the histogram is characteristic of every event which we have 
studied to date. The two left-hand peaks are produced by protons and 
alpha particles respectively, while the riglit-hand peak is due to 
peirticles having ionic charge For C and 0, and 8 5 Q < 12 for iron. 

The apparent energy per charge range for iron is lower because for low 
energy iron the energy measured in a solid-state detector is considerably 
less than its incident energy (Campbell and Lin, 1973). The large dip 
between the alpha peak and the Q = 5-8 peak is statistically consistent 
with a complete absence cf ionic charges 3 and In no event which we 
have examined have we ever seen evidence for the presence of these two 
charge states . 

III. OBSERVATIONS 

Our observations of low energy particles produced in the series of 
flares beginning on 197^ July ** have been reported earlier (Sciambi et al . , 
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197^). We sunaaarize the results of our analysis of that event: (l) Rate 

data for 15-600 keV per nucleon ions heavier than helium indicate / :.rong 
bend-down in the energy spectra below 100 keV per nucleon, even late in the 
event, which cannot be attributed to velocity dispersion; (2) the proton 
and helium spectra show no such bending within the energy range of our 
detectors {i 130 keV for protons and ^ UO keV per nucleon for alphas); 

(3) pulse-height analysis data support the observations listed above, and 
indicate that carbon and oxygen are almost fully stripped, having mean 
charge states of 5*5 and 7.2, respectively; and (U) the average a/CNO ratio 
for the event, as determined from PHA data, is about l6. 

We have investigated a total of eight additional solar particle 

events observed by IMP-8 since its launch in late 1973. These data 

indicate that the features of the CNO spectrum in the July event were not 

unusual. In Figiire 3, for example, we show proton, helimn, and CNO spectra 

ll 

averaged over the time period 197^ June 10 O” to June 12 o". This event, 
associated with a small solar flare on June 8, is typical in that the CNO 
spectrum flattens at low energies. The one unusual feature of this event 
(as compared with the other eight examined) is that the helium also seems 
to show a flattening; in fact, this is the only case in which we have seen 
such an effect in the alphas. 

Figure h shows the time variation of the spectra index y (J "v E ^) 
for heavy ions between 20 and 14^0 keV per nucleon during the 10 June event. 
The effect of velocity dispersion is evident up to the beginning of day 
162, at which time the spectral index settles down to a low, more-or-less 
constant value and remains at that value to the end of the event. The 
bend-over in the spectrum (low spectral index) cannot, therefore, be 
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attributed to simple velocity dispersion. 

Using the arrival times of the low energy heavy ions of various energies, 

it is possible to calculate their diffusion coefficients, ic, based on travel 

distances (along the average spiral magnetic field) of 1.2 AU. We obtain 

20 20 2 

values of < of between it x 10 and 6 x 10 cm per sec, for ions with 

energies between 20 to 250 keV per nucleon. The very small (<<3^) residual 
anisotropies which we observe after correcting for the Compton-Getting 
effect (Ipavich, 197^*) are consistent with this low value of the diffusion 
coefficient. 

Table 2 summarizes the observations for the nine particle events, all 
of which are associated with solar flares in the western limb of the sun. 

Four of these events show an "imusually low” spectral index (y < 1) at low 
energies for CNO, while, with the exception of the June 10 event, none show 
a clear alpha turndown. We must caution, however, that since for alpha 
particles the low energy limit of our observations is 4o keV per nuclton, 
a factor of two higher than the lowest heavy ion energy, the possibility of 
an alpha turndown below 1+0 keV per nucleon cannot be ruled out. On the 
other hand, one can definitely say that if a turndown is present in the 
alpha spectra, it is not nearly as strong as for the case of the heavy ions. 
Elemental abundance information from ULET indicates that for the nine events 
observed, the abundance of nitrogen relative to (carbon and oxygen) was <10/5 
at 500 keV per nucleon for each of the cases examined. It is therefore 
reasonable to assume that very little nitrogen is present at the lower energies 
as well, and we are thus able to make determinations of the ionization states 
of carbon and oxygen. We find that in all event,s the ionization states of 
carbon and oxygen are about 90 % of the nuclear charge, i.e., these ions 
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are nearly fully stripped. The a/CNO ratio at Uo keV per nucleon shows a 
great variability from event to event, hut there is no correlation of CNO 
ah\mdance with the degree of turndown in the spectra; nor is there a 
correlation with helium flux, CNO flux, or location of the flare. There 
is, however, a strong correlation of CNO spectral index Yqjjq vith the 
low energy (s 150 keV) proton flux averaged over the event. This 
interesting correlation is evident from Figure 5, where we plot the CNO 
spectral index from U0-80 keV per charge ^ the "event averaged" 130-220 keV 
proton flux. The trend is clear: the larger the flare as measured by 

the low energy proton flux, the stronger the bend in the spectrum. 

IV. DISCUSSION 

After passage throu^ a sufficient amount of "cold" matter, the 
charge states of ions approach limiting values Q* given by the semi-empirical 
relation (Betz, 1972) 

Q* = Z[1 - 1.032 exp (-137B/Z*^^) ] 

where Z is the atomic number and $c the velocity of the ions. For the 
energies observed by us, Q* is computed to be about 2.5 carbon and 8.9 
for oxygen. These numbers are clearly inconsistent with our results. 

The heavy ions we observe may have been stripped through impact 
ionization by keV electrons produced in the flash phase of the flare, as 
suggested by Braddy et ^ (1973). However, the charge states, while high, 
do not indicate complete stripping, as we would expect if this mechanism 
were the dominant one. 

We note that from the work of Jordan (1969) and others (Allen and 
Dupree, 1967; Lange and Scherb, 1970) the ionization states observed are 
consistent with those expected for ions in thermal equilibrivim with an 
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electron gas at a temperature of 1-2 x 10 ®K. In addition* the charge states 
are very similar to those observed in the solar wind {Bame ^ » 19T0). 

We thus suggest that the low energy heavy ions seen near earth have 
achieved charge state equilibrium in regions of the solar corona from 
which the solar wind originates. 

We now address the question of what causes the observed tiarndown in 
the CNO spectrum. The first possibility which comes to mind is the effect 
of ionization energy loss in xhe spectra during propagation from the 
acceleration region, and/or storage in the corona. In this case it is 
important to use the correct form of the energy loss equation for low 
particle energies. The most commonly used form, dE/dx 'v l/E, breaks down 
as the particle velocity approaches the electron thermal velocity 
For an electron temperature of ^-lO^^K, this corresponds to a particle 
energy of a,225 keV per nucleon. The problem of the slowing down of an 
ion in a fully ionized plasma has been treated by several authors. We 
shall use the results of Itikawa and Aono (1965), who find a form of 
the energy loss equation which is valid for all values of the ion velocity: 


dE. 
1 

dx 


= - 1.512 X 10 


-12 \ ^ 

T M. 
e 1 



3 \l/2 


+ 9.13 


erfM- 1.129 

2 u 

u 


( 1 ) 


with u = T9. 5 /e^/T^, where E^ is the ion energy in keV per nucleon, T^ is 

the electron temperat\ire of the medium in ®K, is the ion mass in Ah-fJ, 

X is the path length in cm, n is the electron density of the medium in 

e 

-3 2 /"^ 2 

cm and erf is the normal error function: erf(u) = — / exp(-x ) dx. 

/ff *^o 

In the low energy limit << the specific energy loss, dE/dx is 

proportional to or while in the high energy limit 


9 


(▼. >> V.. ), dE/dx l/E. . We assume that the particles were initially 

ion the 1 

accelerated to a power law spectrxim in kinetic energy, J(E) ■ J^E and 
ask what effect ionization loss could have cn the spectrum. 

It can he shown (Ipavich and Sciambi, in preparation) that propagation 
of all particles through the same amount of material causes the spectral 
index, for low energies, to approach a limiting value of 1/2, regardless of 
the input spectral index. Since for some events there is evidence that y 
is less than 1/2, we look to a combination of propagation and storage to 
produce the observed effect. 

In Figure 6 we show numerical integrations of equation (l) for various 
ions, for an initial power law spectrum J(E) E”^, propagated outward 
from an injection point at 1.2 R_ to a storage area at 2 R_, stored at 2 R„ 
for 10 hours, and then propagated out to earth. The corona has been assumed 
to be isothermal, with an electron temperature of 1.5 x 10^°K. The particular 
values of the storage radius and time do not affect the limiting slopes of 
the curves in Figure 6, only how quickly these limits are reached. The 
factor Q /M in equation (l) is apparent in the difference in degree of 
flattening for alphas and heavier nuclei. Clearly a ”kink” occurs in the 
spectrum. However, it is at an energy about ten times higher than observed 
(see Figure 3). In addition, the energy spectra of both alphas and protons 
should also show a bend-over, a feature which is not evident in the data. 

The energy at which the bend occurs is, to first order, a f\inction 

only of the paiometer u in equation (1) , or of E, /T . In order to make 

1 e 

the calculated spectra consistent with observation, the bend energy must 
be lowered by about a factor of 10. This can be accomplished by lowering 
the temperature of the storage by a factor of 10. But from the observed 
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char^ states of the ions ve may put limits on the te&qierature of the 
stori^ region. Specifically, from the upper limits of the ahun^mice of 

S 

C observed, the storage temperature must have been greater than 8 x 10 ®K 

(Jordan, 1^9 )t f&P too high to produce the flattening at the observed 

energy. We my also rule out storage at low temperatures (with low 

equilibrium chsu^ge states) followed by stripping for the following reason: 

A turndown in the spectra of CHO which is more pronounced than for alphas, 

2 

requires that in the storage region the parameters in luation (1) must 

£ 

have been greater for CNO than for alphas. Since for alphas, n- * 1, 

2 

Q > 12, or Q > 3.5 for carbon, which i^lies that in order for these hi^ 
charge states to be maintained, the temperature in tl» storage regions 
must again have been > 8 x 10^®K. Therefore, we argue that coronal 
propagation ^d storage alone could not have produced the spectral 
flattening at the energies where we observe it . 

Verzariu and Krimigis (1972) have pointed out that by allowing the 
diffusion coefficient k to depend on particle velocity, a bend in the 
spectrum can be produced by adiabatic deceleration. But, since the 
energy change mechanism would depend only on particle velocity, different 
species observed at the same velocity should have experienced identical 
propagat ion effects, and should have “mapped back" to the same source 
velocity. However, in the present case, the absence of a strong alpha 
turndown at the same kinetic energy per nucleon (or velocity) as for carbon 
and oxygen makes it unlikely that the bending is caused by a velocity- 
dependent n. Using an identical argument, we rule out the effect of a 
rigidity-dependent diffusion coefficient, since by virtue of their high 
state of ionization the carbon and orygen have nearly the same rigidity as 




the alphas at the saute energy per nucleon. Adiabatic deceleration can 
certainly char.ge the shap‘ of the spectra; our point here is that this 
decelez%tion should produce Identical changes in the spectra of all 
species observed. 

We suggest that the observed bend in the spectra is caused In the 
following manner; particles (alphas and heavier nuclei) are accelerated 
in a source region to a power law spectrum In kinetic ener^. Ihe 
particles are then stored in a region of coronal temperatures (T^'v 1.5 x 10^®K), 
which may or may not be the same as the source region. Ibis trapping region 
is probably within 3 of the solar surface (Anderson, 19T2). Here the 
particles lose energy through ionization loss, their spectrum flattens 
at about 800 keV per nucleon and some of the ions are stripped of their 
outer electrons. From our observation of the observed correlation between 
the degree of flattening and the low energy proton flux, we conclude that 
the storage time is proportional to the "size” (i.e., total energy) of the 
flare. The particles are then diffusively propagated outward to the 
earth, losing energy by adiabatic deceleration. Here we assume that the 
diffusion coefficient is independent of particle parameters, so that the 
adiabatic deceler.'i'ion does not affect the shape of che spectrum, but 
simply iwDVes it to lower energies (Venkatarangan and Lanzerotti, 1975 ;• 

Under these conditions, the required energy losses may be produced by a 
diffusion coefficient cf the order of 5 x 10 cm sec (Parker, 1965 ) which 
is consistent with our observations. Thus, th*» "bend” in the spectrum of 
CNO is BK>ved down to the energy at which it is observed. Dtatistical errors 
in our helium observations prevent an unambitious identification of the 
predicted bending in the helium spectra. 
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Because theoretical results derived for alpha particles apply 
equally for protons, there should also be a bend in the proton spectrum 
to match that of the alphas. Since the lowest observed energy for protons 
is even higher, than that for alphas, we do not expect to see this bend. 

We also note that all species appearing in Figure 6 were initially equally 
abundant. After storage, the heavy ion abundances are depressed relative 
to helium. Thus, in order to explain events which show enhanced heavy 
element fluxes, it is necessary either to have . ’’super enriched” source, 
or to invoke some form of preferential acceleration. 

Althou^ we believe that we have proposed a reasonable mechanism for 
producing most of the observed features of the carbon and oxygen spectra, 
our conclusions at this point are tentative. Additional event statistics 
may shed further light on the relative contributions of coronal storage 
and energy losses in producing the observed charge states and spectra. 

V. SlMlARi 

In a total ol nine solar particle events of various sizes we have 
observed a bend-down of the energy spectra of carbon and oxygen below 
100 keV per nucleon. In addition we find that: 

(a) In all events examined carbon and oxygen nuclei are nearly fully 
stripped; the observed ionization states are consistent with 
thermal equilibrium of these ions with the medium at a temoerature 
of a.lO^°K. 

(b) There is a great variability from event to event In the a/CNO 
ratio at Uo keV per nucleon. 

(c) There is, in all but one case, no clear helium bend-down. 
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(d) There is no apparent correlation of the degree of turndown 

in the CNC spectrum with the relative CNO abundance, the helium 
flux, the CNO flux, or the location of the flare. 

(e) The amount of bending In the CNO spectrum is strongly correlated 
with the event-averaged flux of protons in the energy range 


130-220 keV. 

These observations are consistent with storage of an initial power 
law spectrum at coronal densities of '^3 x 10 cm”^ and temperatures of 
1-2 X 10^°K for 10 hours, followed by strong adiabatic deceleration. 
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Table I 

Energy Ranges and Charge Responses for IMP-8 Detectors 


Detector 

PHA Charge Threshold 

Energy Range 
(keV per charge) 

Charge 

Thresholds 

Sectored 

PI 

Q>6 

37-55 

Q> 6 

Yes 

P2 

Qt.2 

65-100 

q>2/Q>6 

Yes/Wo 

P3 

QU 

130-220 

Q> 1 

Yes 

P4 

Qil 

160-230 

Q>1 

No 

P5 

background 

>160 keV 

— 

Yes 

P6 


390-600 

q>i/q^6 

Yes/Yes 

P7 

Qil 

640-1100 

Q>l/Q>2 

Yes/No 

P8 

QRl 

740-1200 

Q>1 

Yes 

N 

electrons 

600-800 keV 

— 

Yes 


Characteristic of Nine Solar Flare Particle Events 
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Figure Captions 

Figure 1 — Cross-sectional view of the University of Maryland 
electrostatic deflection sensor as flown on the IMP-6 spctcecraft. The 13 
solid state detectors are surrounded by a plastic scintillator anti- 
coincidence. Six of the detectors (open rectangles) are used in an anti- 
coincldene® mode to prevent analysis of penetrating particles. 

Figure 2 — lypical pulse hei^t histogram from the IMP-8 experiment, 
obtained from detector P3 during a solar particle event. 

a) shows the histogram before correction for background (the smoothed 
background curve is shown as a dashed line). 

b) show the histogram after subtraction of backgroimd. Note the three- 
peaked form of the corrected histogram. 

Figure 3 — Proton, helium, and CNO group energy spectra in the 
10 June lirT^* event. 

Figure U — CNO spectral index (h0-80 keV per charge) ^ time for the 
10 Jvine 197^ event. 

Figure 5 — CNO spectral index (U0-80 keV per charge) v£ "event- 
averaged” proton flux (130-220 kev) for the nine particle events observed. 
Note the strong correlation with event size. 

Figure 6 — Calc\ilated in differential spectra after coronal storage. 
The curves were obtained by numerically integrating equation (l) outward 
from the injection point. The strai^t line is the injection spectrum 
for all species. 
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